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Contributions of two-pseudoscalar and vector-pseudoscalar meson loops to the p-uj mass splitting 
are evaluated in a covariant model based on studies of the Schwinger-Dyson equations of QCD. The 
role and importance of the different time orderings of the meson loops is analyzed and compared 
with those obtained within time-ordered perturbation theory. It is shown that each meson loop 
contributes less than 10% of the bare mass, and decreases as the masses of the intermediate mesons 
increase beyond approximately mp/2. A mass splitting of m„ — rup « 25 MeV is obtained from the 
nn, KK, utv, pn, cjrj, prj and K*K channels. The model is then used to determine the effect of the 
two-pion loop on the p-meson electromagnetic form factor. It is shown that the inclusion of pion 
loops increases the p-meson charge radius by 10%. 



I. INTRODUCTION 

In a fully-interacting theory, the number of con- 
stituents of a hadron is not conserved, and the picture 
of a hadron as being comprised of three quarks (baryon) 
or a quark and antiquark (meson) is naive. One way to 
go beyond a picture in which hadrons are simple bound 
states of valence quarks is to consider their mixing with 
multiple-hadron states. The success of models which de- 
scribe hadrons in terms of only a few degrees of freedom, 
such as the constituent-quark model, suggests that mix- 
ing with multiple-hadron states should represent only a 
small correction to the predominant valence quark struc- 
ture of the hadron. 

Some observables may be sensitive to contributions 
that arise from mixings of hadrons with multiple-hadron 
states. These observables would provide a means to 
probe the structure of hadrons beyond that of the va- 
lence quarks. One such observable is the mass splitting 
between the p and lu mesons, which is the subject of this 
article. 

The and uj have a similar quark substructure, dif- 
fering only in their total isospin (and hence G-parity). 
SU(3)-flavor violating contributions due to the different 
masses of the u, d and s quarks tend to split the masses 
of the p and to mesons. However, these contributions 
alone do not account for most of the mass splitting [Q. 
Other contributions arise from the mixing of the p and lu 
mesons with multiple-hadron states of differing G-parity. 
An important example of this is the shift in the mass of 
the p meson due to its mixing with the two-pion state. 
As mixing between the u> meson and the two-pion state is 
forbidden by G-parity, the mass of the to meson doesn't 
receive such a contribution from the two-pion state. The 
difference between the p and oj mesons due to mixing 
with all possible multiple-hadron states results in an ob- 



servable splitting of the masses nip and m^j, which would 
otherwise be degenerate (in the absence of SU(3)-flavor 
breaking). 

In Ref. 1^, the mass shift of the p meson due to its 
mixing with the two-pion state was investigated using 
an effective chiral Lagrangian approach. A form factor 
was introduced into the resulting dispersion relations for 
the p-meson self energy n^,y(g^) to render it finite. Sev- 
eral different forms of this form factor were considered 
and each produced similar results. For example, using 
a dipole form factor, it was found, over a large range of 
dipole widths 0.6 < A^ < 4.0 GeV^, that the p-meson 
mass is shifted by roughly —10 to —20 MeV due to its 
coupling with the two-pion state. 

The purpose of the study of Ref. |^ was to determine 
whether results obtained by lattice studies of QCD for the 
p meson mass are significantly affected when the p meson 
is "unquenched" ; that is, when the p meson is allowed to 
mix with two-pion states. Such mixings with multiple- 
hadron states are typically neglected in lattice studies. It 
is clear from the study, which relies on very little dynam- 
ical modeling, that the inclusion of intermediate meson 
states into the vector-meson self energy Ilp^^{q'^) is indeed 
a small effect which can be treated as a correction to the 
predominant quark-antiquark structure of the vector me- 
son. 

An extensive study of the p-tu mass splitting is carried 
out in Ref. In this work the contributions of many 
different two-meson intermediate states to the masses of 
the p and u are evaluated within a nonrelativistic frame- 
work, using time-ordered perturbation theory to evaluate 
the loops. The three- meson vertices and their momentum 
dependence are evaluated using a string-breaking picture 
of the strong decays based on the flux-tube model of the 
confining interaction, and quark pair creation with vac- 
uum i^Pn) quantum numbers. It is noted that the p-ui 
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splitting is directly proportional to an OZI-violating mix- 
ing between ground-state vector uu and dd mesons. This 
mixing is due to two-meson intermediate states which can 
link the initial and final mesons. It is shown in Ref. 
that, in the closure limit where the energy denominators 
for the intermediate states are constant, the ^Pq pair- 
creation operator cannot link these two vector states, so 
that the OZI rule is exact in this limit. The presence of 
a small OZI-violating mass splitting is, therefore, inter- 
preted as a deviation of the sum over intermediate states 
from this closure limit due to the differing properties of 
the physical intermediate states. 

Although the authors of Ref. point out that their 
calculation of this mixing is not formally complete due 
to the absence of pure annihilation diagrams where the 
intermediate state contains only glue (which can be writ- 
ten as a sum over heavier glueball states and are likely 
suppressed), the sum over intermediate states is much 
more extensive than that carried out here. The conver- 
gence of this sum, although convincingly demonstrated, 
is rather slow in their calculation, requiring many dif- 
ferent meson pairs with many possible relative angular 
momenta. This can be illustrated by examining the size 
of individual terms in the sum; for example, the tttt loop 
contributes —142 MeV to the p mass, and the utir loop 
contributes —146 MeV (the corresponding pir loop con- 
tributes roughly three times this amount to the uj mass). 
It is pointed out in Ref. |^ that without a flux-tube over- 
lap function and a form factor at the pair creation ver- 
tex, both of which are necessary for a theoretically con- 
sistent picture of the strong decays and suppress high- 
momentum created pairs, this sum would not formally 
converge. Note that neither of these features of the strong 
decay model are essential to fitting the strong decays of 
mesons in this picture, but were essential to the conver- 
gence of the loop expansion in Ref. A global fit to 
the meson decays was carried out in order to fit the de- 
cay model (in particular the pair creation vertex form 
factor). 

It was concluded in Ref. |^ that reasonable results 
for the p-LU mass splitting could be obtained only when 
the effects of many two-meson intermediate states are in- 
cluded in the loop expansion. Such a conclusion raises 
questions about the validity of the usual picture of a me- 
son as being well described as a bound state of a valence 
quark and antiquark. 

The motivation for carrying out the present study is 
to attempt to reconcile the small effects on the p mass 
due to the two-pion loop found in Ref. |^ with the more 
substantial mass shift found in Ref. and to examine 
the consequences for the convergence of the sum over in- 
termediate states performed in Ref. [||. The p-oj mass 
splitting is evaluated using meson transition form fac- 
tors obtained from a covariant, quantum field theoretic 
model based on the Schwinger-Dyson equations of QCD. 
A similar approach was employed in Ref. |0] for the tttt in- 
termediate state. The present study extends the work of 
Refs. by considering the additional channels KK, 



LUTT, pTT, ijjr], prj, and K*K. (In Ref. ||], a non-local 
quark four-point interaction was used to study the p- 
LO mass splitting and it was shown that the three-meson 
intermediate state tttttt produces a negligible mass shift. 
Hence, three-meson intermediate states are neglected in 
the present study.) The model employed is fit to vari- 
ous properties and electromagnetic form factors of the tt 
and and has been shown to give reasonable results for 
the decays and electromagnetic form factors of the light 
vector mesons in Ref. . 

Although a very different framework is employed here, 
results that are qualitatively and quantitatively similar 
to those found in Refs. are obtained for the p- meson 
mass shift due to the two-pion loop. In each of the 
seven channels considered here, a mass shift of less than 
10% of the total mass of the p is observed. These mass 
shifts partially cancel each other when calculating the 
difference rrit^ — nip, as observed in Ref. ||^, and we ob- 
tain a net mass splitting of — nip « 25 MeV. The 
experimentally observed value for this mass splitting is 
- mp = 12 ± I MeV. 

Of course, the calculation of the p-u) mass splitting de- 
scribed herein is incomplete. Other channels, such as the 
two- vector-meson intermediate states, might be expected 
to contribute significantly H). Nonetheless, these results 
are consistent with those obtained in Ref. and in- 
dicate that the sum over two-meson intermediate states 
should converge quite rapidly as the masses of the in- 
termediate state mesons increase, in contrast to what is 
found in Ref. 

The dependence of the vector meson self energy on the 
masses of the intermediate mesons is also explored in de- 
tail. It is shown that contributions from two-meson loops 
with meson masses above nip/ 2 contribute only small 
amounts to the mass shift, and their contributions de- 
crease rapidly as the masses of the intermediate mesons 
are increased. This is the essential mechanism for the 
expected rapid convergence of the sum over two-meson 
loops and is insensitive to the details of the model meson 
transition form factors, the only dynamical input em- 
ployed in the present calculations. 

Another observable that is sensitive to the inclusion of 
meson loops is also explored in this study for the first 
time. With the same p tttt transition form factors 
used in the calculation of the p-uj mass splitting, the 
contribution of tttt loops to the p-meson electromagnetic 
form factor is calculated. This model predicts a 10% 
increase in the charge radius of the p meson relative to 
that obtained in Ref. ^ which considered only the quark- 
antiquark structure of the vector meson. The inclusion 
of pion loops leads to a similar 10-15% increase for the 
TT- meson charge radius The results from the present 
study also indicate that of the two-meson states consid- 
ered, the two-pion state provides the most significant cor- 
rection to the vector meson charge radius. 

The brief lifetime of the p meson makes a direct mea- 
surement of its electromagnetic charge radius impossible 
at present. However, the size of the neutral vector me- 
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son can be extracted from high-energy data for diffrac- 
tive, vector meson photoproduction on a nucleon. It is 
possible, usinga semi-empirical model of diffractive pho- 
toproduction 10] , to relate the diffractive radii of mesons 
to their charge radii. Given the dominance of the two- 
pion state this implies that the increase in the charge 
radius of the p meson due to the two-pion loop corre- 
sponds to an analogous increase of the diffractive radius 
of the p meson. Since G-parity forbids the uj meson from 
receiving such contributions from the two-pion state, a 
measured difference in the diffractive radii of the p and 
LD mesons can provide a means to probe the magnitude 
of the two-pion contribution to p meson observables. 

T he o rganization of this article is as follows. In 
Sec. II A , the general Lorcntz and flavor structure of the 
meson transition amplitudes, necessary for the calcula- 
tion of the vector-meson self energy and electromagnetic 
form factors, are derived. These elements are then em- 
ployed in a study of the Schwinger-Dyson equation for 
the vector meson propagator in Sec. [IB. The mass of 



the vector meson is related to the real part of the vector- 
meson self energy, and the total decay width of the vec- 
tor meson is related to the imaginary part. These same 
elements are then employed in a study of the p meson 



electromagnetic form factors in Sec. [IC. 

In Sec. m, the contributions to the vector-meson self 
energy arising for the different time orderings of the in- 
termediate meson propagators are investigated. Aspects 
of time-ordered quantum mechanical frameworks, such 
as that employed in Ref. and the Lorentz-covariant 
Euclidean-based quantum field theoretic framework em- 
ployed here are discussed and contrasted. 

In Sec. IV, the dynamical model used to calculate 



the off-mass-shell behavior of the meson transition am- 
plitudes is described. The model describes the meson 
transitions in terms of loops involving nonperturbatively- 
dressed quarks and model Bethe-Salpeter amplitudes 
that describe the quark-antiquark substructure of the 
mesons. The dressed-quark propagators and Bethe- 
Salpeter amplitudes employed herein were developed and 
tested in numerous studies of meson observables |^,||j9| 
and in numerical studies of the quark Schwinger-Dyson 
equation of QCD |o|-0. 

Result s for the p- an d (jj-mesonself energies are given 
in Sees. VA and VB. In Sec. VC, numerical results 
for each of the possible time orderings of the two-pion- 
loop contribution to the self energy are provided. Results 



study are given in Sec. VI 
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II. MESON-LOOP CONTRIBUTIONS 

In this section some of the formalism necessary to de- 
scribe the effects of two-meson loops on the vector-meson 
self energy and electromagnetic form factors is presented. 



In Sec. [I A, the most general Lorentz-, parity- and 
SU(3)-flavor-covariant amplitudes that describe the cou- 
pling of a vector meson to two pseudoscalar mesons or 
a vector and a pseudoscalar meson are presented. The 
transition amplitudes are written in terms of a coupling 
constant and an off-mass- shell transition form factor. 
Both are defined so that the transition form factor is 
equal to one when all three mesons are on shell. In the 
case of the p meson coupling to two pions, the coupling 
constant associated with the transition amplitude is re- 
lated to the decay p tttt and is found by experiment 
to be gpTTTT — 6.03. The other coupling constants consid- 
ered herein as well as the off-mass-shell transition form 
factors are obtained from considerations of SU(3)-flavor 



invariance or the dynamical model described in Sec. IV 



In Sec. [IB, the Schwinger-Dyson equation for the vec- 
tor meson propagator is given in terms of the vector me- 
son self energy Ilfj,^{q). Expressions are derived that give 
the self energy in terms of loop integrations involving 
the propagators for the two intermediat e me sons and the 
transition amplitudes described in Sec. [I A. The imagi- 
nary part of the vector-meson self energy is related to the 
total decay width of the vector meson, and the real part 
of the self energy provides a contribution to the vector 
meson mass. 



In Sec. |II Cj the contribution of the two-pion loop to 
the p-meson electromagnetic form factors is considered. 
A correction to the electromagnetic (EM) charge radius 
of the p meson due to pion loops is obtained, and is com- 
pared to the results obtained from a study of the quark- 
antiquark contribution from Ref. Q. It is argued that 
since G-parity forbids the lo meson coupling to two pi- 
ons, and since the two-pion contribution dominates, mea- 
surement of the difference of the p meson and co meson 
charge radii provides a means to determine the effects 
of including such intermediate meson loops. The elec- 
tromagnetic form factors are given in terms of the same 
meson transition amplitudes used in Sec. 11 B to calculate 



the self energies and decay width of the p meson. Hence, 
the calculation of the contribution of pion loops to the 
p meson EM form factor provides a further test of the 
self-consistency of the present approach. 



A. Meson transition amplitudes 

The most general coupling of a p meson to two pseu- 
doscalar mesons (tt or K) can be written in terms of the 
following action: 



d^xd\d'z7r\x)7r^{y)pl{z)A^^\x-z,y~z), (2.1) 



where p^^{z) is the p meson field with flavor index k and 
Lorentz index p, and 7r*(a;) is the pseudoscalar meson 
field with flavor i. The nonlocal three-point coupling is 
described by the amplitude 
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which can be written in terms of a coupUng constant 
gijk, where i, j, and k are SU(3)-flavor indices, and 
a form factor f^^^ {pi,P2) which is a function of the 
Lorentz invariants pf, p\ and = {pi + ^2)^- The 
form factor f^^^{pi,P2) is defined to be equal to unity, 
when aU three mesons are on their mass shell; that 
is, f^^^{pi,P2) = 1 when pj ^ p^ ^ -ml, = 
{pi + P2)^ = ^iTT-'i- Hence, one may write: 
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Arbi>P2) = ^(Pi-P2)p g''" r^^bi,P2). (2.3) 



In general, Eq. ( p.3| ) may also have a term which is pro- 
portional to {pi +P2)p,- Even if such a term were present, 
it could not contribute to the self energy of an on-shell 
vector meson since its contraction with a spin-1 polariza- 



tion vector is zero. For p tttt. 



ijk 



9p^ 



ijk 



From the above amplitude, one obtains the invariant 
Feynman amplitude for the decay of, say, the into two 
charged pions: 



+ {p,)-^-{p2)\T\p\q,\)) 

SffpTTTT Pi A), 



(2.4) 
(2.5) 



where Sfi{q, A) is the polarization vector for a p meson of 
momentum q and helicity A. The resulting p — > tttt decay 
width is 



47r 12 
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4to; 



3/2 



(2.6) 



The meson transition form factor f^^^ {pi,P2) does not 
appear in the relation for the decay width because it is 
unity when all mesons are on their mass shell. Its value 
away from this point is therefore not directly observable 
and must be calculated from a dynamical model. 

Similarly, the coupling of an uj meson to a p and tt 
meson is described by the action: 



d^x d% d^z pI{x)t:' {y)uj^{z)K%Xx~ z,y- z). (2.7) 



Here, uju{z) is the w-meson field and A^^(a; — z, y — z) 
is the nonlocal transition amplitude describing the cou- 
pling of two vector mesons to a pseudoscalar meson. In- 
troducing the coupling constant g^pTr and form factor 
f^^^{pi,P2), considerations of Lorentz covariance give 
the Fourier transform of A^^(a; — z, ?/ — z) as 

^Up^,P2) = e^ucp^^^ 5''gu>p.f''''{puP2) (2.8) 



where pi and p2 are the momenta of the p and tt mesons, 
respectively, and f^^^{pi,P2) = 1 when q^ — {pi + 



and P2 
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-to: 



Unlike the 



P2r = -mZ, Pi = -to;, 
previous case, a lack of phase space prevents the decay 
Lu — !■ pTT, SO that it is impossible to determine the value 
of the coupling constant gcop-ir from experiment. Hence, 
both the meson-transition form factor f^^^ {pi,P2) and 
the coupling constant gujp-jT must be determined from a 
model calculation. 

In Sec. IV, a dynamical model is described and used 
to determine the transition form factors f^^^ {pi,P2) 
and f^^^ {pi,P2), and the coupling constant g^jp^r- The 
model has been shown to give good results for the simi- 
lar, physically accessible decays p ^ ttj and K* K'y, 
in which a vector meson decays into a vector and pseu- 
doscalar 1^. The value of the coupling constant gpTm- 
— 6.03 used herein is extracted from the e xper imentally 
observable p —^ tttt decay width using Eq. (2.6). 



B. Vector-meson self energy 

The spin-1 nature of a vector meson entails the follow- 
ing form for the vector-meson self energy when the vector 
meson is on mass shell (g^ — —my): 



H,.(g)=(V-^)n(9^) 



(2.9) 



Away from the on-mass-shell point q^ = ~~rny, the 
vector-meson self energy is of the general form: 



Tlp.iq) 



q^ J q 



n(,2) + M!ie(,2). (2.10) 



To project out the Lorentz scalar function Tl{q^), a pro- 
jection operator ^T^i,{q,q^) where 



T^_iu{q, -TO^) = + qpqv/m^ 



(2.11) 



is employed. 

In the following, the dressed vector meson propa- 
gator D^v{q) = Tf^^{q,-my)b{q^) is given in terms 
of the "bare" vector meson propagator Df^^{q) = 
Tp,y(g, — TOy)Z?o('Z^) and the vector-meson self energy 
Hp,y(g) — Tf^„{q,—my)U{q'^) which describes the mix- 
ing of the vector meson with meson loops (and in general, 
multiple- hadron intermediate states) . The "bare" vector- 
meson propagator DQ^q^) is presumed to arise from the 
quark-antiquark structure of the vector meson, and does 
not include the effects of multiple-meson intermediate 
states. Defined in this manner, quark confinement re- 
quires that the "bare" vector-meson propagator Do{q^) 
is a real- valued function of q^. The presence of a non- 
zero imaginary part for DQ{q^) would entail a loss of flux 
for vector-meson propagation, which is a signal for the 
existence of open decay channels for the vector meson. 
In our approach, such decay channels arise from the cou- 
pling of the vector meson with multiple- hadron states. 
Mixing the vector meson with multiple-hadron states re- 
sults in a dressed, vector-meson propagator D{q^) that is 
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a complex-valued function of the vector-meson momen- 
tum . 

The dressed, vector- meson propagator D{q^) is ob- 
tained by solving the Schwinger-Dyson equation: 



D = Dn 



DqUD. 



(2.12) 



This integral equation gives the dressed vector meson 
propagator D{q^) in terms of the bare propagator Do{q^) 
by dressing the vector meson by two-meson loops as given 
by n(g^). As an integral equation, Eq. (2.12) necessarily 
implies the insertion of infinitely many self-energy n(q^) 
contributions into the dressed, vector-meson propagator 
D{q^)': hence, the dress ing o f the vector-meson propaga- 
tor resulting from Eq. ( 2.12 ) is inherently nonperturba- 
tive. 

The solution to the Schwinger-Dyson equation is 



D = Do + DallDa + DanDoUDo 
1 



Dn 



I + UDo 



so that 



D-^ = (1 + UDo)D^^ = D^^ + n. 



(2.13) 



(2.14) 



Near q^ = —my, the difference between the dressed 
and bare inverse propagators is given by 

Uiq') = D-\q')-D^\q') 

= (1 + SZ) {q^ + m\ - imyV) - {q^ + m^), (2.15) 

where TOq is the bare mass of the vector meson, F is the 
total decay width of the vector meson, and 1 -I- 5Z is 
the field renormalization for the vector meson, chosen to 
ensure that the asymptotic vector meson has unit nor- 
malization. In Sec. [V, the Bethe-Salpeter amplitudes. 



which describe the quark-antiquark substructure of the 
meson, are normalized to unity. Thus, from the defini- 
tion of the normalization 1 -I- 5Z , as given by Eq. ( 2.15 ), 
it follows that the correction to the normalization that 
results from the inclusion of meson-loops is given by 5Z . 
For meson loops to be considered as a small effect, 5Z 
must be smaller than one. 

Comparison of Eq. ( p.l5[ ) and Eq. (2.14) gives 



Im(n) = (1 + 5Z) mv F 
Re(n) = 5Zq^ + (1 + 5Z)m 



(2.16) 
(2.17) 



Eq. (2.17) provides a means to determine the normaliza- 
tion 6Z from the self energy Ii{q'^): 



from which it follows: 



= ml+ Ren(- 



= 5Z, 



4)- 



(2.18) 



(2.19) 



This equation gives the mass of the vector meson my, 
in terms of the bare mass mg (arising from the quark- 
antiquark substructure of the vector meson) and the self 
energy Il{q^). It can be solved iteratively, 

my ^ ml + Ren {-ml ~ Ren(-TOy)) 

w + Rcn(-m^) - Ren'(-mg)Ren(-m^), (2.20) 



where the expansion assumes small my — m^. This re- 
lation yields a first-order correction to the vector meson 
mass: 



Ren(-m^). 



(2.21) 



From Eq. (2.18|), one identifies Ren'(— my) with 6Z. 



Thus, the second order correction to the vector meson 
mass shift is — 5ZRen(— m§), which is suppresse d by a 
factor 6Z compared to the first-order term in Eq. ( ^.20 ). 
The importance of higher-order corrections can be esti- 
mated from the momentum dependence of n((7^) in the 
vicinity of TOq- is found to be sufficient (se e Sec . 0) to 
include only the lowest-order term from Eq. ( 2.20| ) to ob- 
tain a good estimate for the dressed, vector-meson mass 
my- 

As stated above, the self energy Il{q^) describes the 
dressing of the vector meson due to two-meson loops. It 
can be described in terms of a fluctuation of the vector 
meson into a multiple-meson intermediate state. This is 
given in terms of meson transition amplitudes, the form 
of which was described in Sec. II A, and propagators for 
the intermediate mesons. The contribution to the vector- 
meson self energy Tl{q'^) due to a two-pseudoscalar- meson 
loop is then given by 

^lUl) =-2 I^Af (-pi,^2)A(pi,mp) 
xA{p2.mp,)Al^'''{pi,P2) 



(2.22) 



where mp and mpi are the masses of the pseudoscalar 
mesons, Pi — \q + k, P2 = \q — k, and the meson transi- 
tion amplitude A^'^(pi,p2), given by Eq. ( |2.3| ), describes 
the transition V PP from a vector meson V to two- 
pseudoscalar mesons P and P with momenta pi and p2, 
respectively. The propagation of the intermediate pseu- 
doscalar mesons is given in terms of the following scalar- 
boson propagator: 



A{pi,mp) = {pI + mp - ie)' 



(2.23) 



where e is an infinitesimal positive number included to 
ensure the proper boundary conditions for the decay 
V — > PP. This process is depicted in Fig. |l], in which the 
meson transition amplitude AJ^'^(pi,p2) is give n in terms 
of a quark-loop triangle diagram. In Sec. |V|, the tran- 
sition form factor is calculated using the dressed-quark 
propagators and Bethe-Salpeter amplitudes of a quantum 
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field theoretic model based on studies of the Schwinger- 
Dyson equations of QCD. In impulse approximation, the 
leading contribution to the meson transition amplitude 
A^'^(pi,P2) is given in terms of the triangle diagram de- 
picted in Fig. |l|. 

Ap plyin g the projection operator ^Tf^^{q,q'^) of 
Eq. ( |2.11| ) to n^^(g) in Eq. I \2.2'^ ) yields the Lorentz- 
scalar vector meson self energy due to two-pseudoscalar- 
meson loops 

{q^ ,mp,mp,) = 



d*k 



A(pi, mp)A(p2, mp>)T^''{k, q), (2.24) 



with 



T''''{k,q)^[k^^ 



k-q^ 



9p..f''''iPi,P2)] . (2.25) 



Similarly, the contribution of a vector-pseudoscalar- 
meson loop to the vector-meson self energy is given by 



d'^k 



rAua(-Pl>-P2)i5^r(pi, my) 



(27r)4 



(2.26) 



where my and mp are the masses of the intermediate vec- 
tor and pseudoscalar mesons, respectively, and the me- 
son transition amplitude A^^(pi,p2) is given by Eq. (2.8). 
The propagator for the intermediate vector meson may 
be written as 



D(,„{p, my) 



A{p,mv) 



Subs titu ting this an d th e relation for A^^j,(pi,p2) from 
Eq. (2^) into Eq. ( 2.26 ), and applying the projection 
operator ^T^„{q,q'^) to the result, yields 



n^^(<Z^, TOy, mp) 

2 r d'^k 
3 



A(pi,my)A(p2,mp).F^^(fc,g), (2.28) 



with 



^^^(fc,g) = 



(fc • qf - k^q^ 



[9ujp^f^^'^{Pl,P2)]'^- 

(2.29) 



Th e flav or structure of the amplitudes in Eqs. ( 2.22| ) 
and (2.26) are given for the particular case of the p 
meson self en ergy. It is easy to generalize the rela- 
tions in Eqs. (2.24) and ( 2.28| ) to obtain the isoscalar 
w-meson self energy. Assuming that the coupling con- 
stants for V PP and V VP are SU(3)-flavor sym- 
metric, one can write the self energy of the p meson 
and the u meson in terms of the Lorentz-scalar func- 
tions U^^{q'^,mp,mpi) and n^^(q^, my, mp) given by 




Pi=k+q/2 




p2=-k+q/2 

FIG. 1. Feynman diagram depicting the quark substruc- 
ture of the vector-two-pseudoscalar meson transition ampli- 
tude A^^'^(pi,p2) (shown as a triangle-shaped quark-loop di- 
agram) given by Eq. (4.1) and the intermediate pion prop- 
agators (dashed curves) necessary for the calculation of the 
vector- meson self energy 11^^ {q^ , mp, mp). 



Eqs. ( 2.24 ) and ( 2.28|) , respectively, for several possible 
channels involving the SU(3)-octet of vector and pseu- 
doscalar mesons. For simplicity, only those two-meson 
loops with the lowest mass intermediate mesons are in- 
cluded in the calculation of the vector-meson self energy 
n((7^). In Sec. ^, it is shown that meson loops involv- 
ing mesons with higher masses are suppressed relative 
to those with lower masses and from Ref. contribu- 
tions from three-meson loops are expected to be very 
small. For the present study, the following channels are 
included: tttt, KK, pir, lot:, prj, uj-q, and K*K. Hence, 
the self energy of the p-meson is given by 



n((j2) = I{PP{q^,m^,m^) + ^11^^ {q^ , jtik , ttik) 

r) + Il^^{q^,mK',rnK) 
,), (2.30) 



^^^(9^m. 



^^^(g^m,p 



(2.27) and that of the oj meson is given by 



n(<7^) -0+i^^^(q^mK,mK) 



-f-3n^^(q^,mp,m„) +Il^^{q^,mK',mK) 



VP I 



\l^P{q' 



(2.31) 



Once the meson transition form factors have been deter- 
mined from a dynamical model, the con tribu tions of the 
self energies can be obtained from Eqs. ( 2.30 ) and ( 2.31 ) 
in a straight-forward manner. 

As an example, consider the calculation of the loop 
integral necessary to obtain li^ ^ (q'^ \ whic h is i llustrated 
in Fig. |l| and is given by Eqs. ( |2.24| ) and (HH), 



nPP{q\mp,mp,) = -^-j 



_(£k_ 
(2^ 



:F^^ik,q) 



[{k + g/2)2 + ml- ie] [{k - q/2)^ + m^, 
Introducing a Feynman parametrization 
1 dx 



aia2 



a2)x + a2\ 



(2.32) 



(2.33) 



the denominators of the propagators in Eq. (2.32) can be 
combined. With a change in variables k k + aq, where 
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a = X ~ 1/2, two of the integrations can be performed 
trivially, giving 



n (q ,mp,mp') 



da 



T {k^aq, q) 



(27r)4 [fc2_(a2_l)q2+^2+^ 



4 1 
"3 (27r)3 



2 / 



T {k — aq, q) 



[B-{a^~\)q'^ + rn? + c 



(2.34) 



where z is the cosine of the angle between k and g, defined 



by fc • g = \/fc^-\ 



g2 z. 



ip/) and m-^ 



{rrip — rrip, ) . An inspection of Eq. ( 2.34| ) reveals that the 
denominator exhibits a zero in the limit that e ^ for 
k^ equal to 



kl = {a^-\)q'~ 



aiv? . 



(2.35) 



This singular behavior leads to a simple pole in the 
vector-meson self energy that is associated with the decay 
of the vector meson int o tw o pseudoscalars. 

The integral in Eq. ( 2.34 ) is evaluated numerically by 
first eliminating one factor of the singularity using 



1 



^2 _ ^2 



1 



k^ ' 



(2.36) 



integrating k^ by parts twice, and noting that Cauchy's 
theorem entails 



lim 



1 



"5+ fc2 - /j2 _ 



fc2 — /Cq 



in5{k'^ 



(2.37) 



where V denotes the princip al par t of the integration. 

The imaginary part of Eq. ( 2.34 ) is obtained quite eas- 
ily. To evaluate the principal part, the fc^-integration do- 
main is divided into two regions. The first region is from 
the origin (fc^ = 0) to twice the distance from the pole 
to the origin (2 k^). The second region is from fc^ = 2 fcg 
to /c^ — Xhe integrations are numerically calculated 
using the method of quadratures. In the first region of 
the integration over fc^, each mesh point of the integra- 
tion is paired with one that is equidistant from the pole 
at fcp to ensure the best cancellation of round-off errors. 
In the s econ d region, the integration is straight-forward. 
In Sec. II A, it is shown that the meson transition am- 
plitudes /^'^^^(pi,P2) decrease rapidly with increasing 
k'^, providing sufficient damping to the integrations in 
Eqs. (2.24) and (2.28) to ensure their convergence. Simi- 



lar techniques are employed in the calculation of the loop 
integral {q"^ ,mv ,mp). 



C. Vector meson electromagnetic form factor 

The emission or absorption of a photon by an on-mass- 
shell p meson is described by the action: 



S = J d^xd%d^z f^^{x) p'p{y) A^{z) 

X e''^ A^^p{x~ z,y-z). (2.38) 

Charge conjugation and Lorentz covariance constrain the 
form of the amplitude Af^apix— z,y— z). Its Fourier 
transform is written in terms of three Lorentz-invariant 
form factors Gi{Q'^) for i = 1, 2, 3, which depend only on 
the square of the photon momentum = {—q — q'Y , 



KMq,q') ^J2^^{Q') Tj:l^{q,q'), (2.39) 



where 



Tj^Uq^q') = {q^-q'^)Ta^{q,-nil)T^p{q',-ml) , 
Tj^Uq^q') = Tf,^{q,-ml)Tps{q',~ml)qs 

-T^piq', -fj4)T„s{q, -ml)q's , (2.40) 



T'Zp{q,q') 



Ta-y(g, -my)q' T: ps{q' , -my)qs , 



my is the mass of the vector meson, and 

rr f 2 \ r , qaqi3 

Tai3[q, -my) = dai3 H —■ 



(2.41) 



One can construct a set of tensors o[*l,o which can be 
used to project out from Af^afj{q,q') each of the form 



factors in Eq. (^; i.e., G,(Q') = 0|:^,3A/.a/3(g, ?')■ 

Alternatively, one may con struc t from linear combina- 
tions of the tensors in Eqs. (2.40) a set of tensors that 
transform irreducibly under spatial rotations in a frame 
where there is no energy transferred to the final p meson, 
so — ((5,0). The form factors associated with this 
set of tensors are the electric monopole Ge{Q'^), mag- 
netic dipole Gm{Q^), and electric quadrupole Gq{Q'^) 
form factors. They are related to the original set of form 
factors by 



Ge{Q') ^ G^iQ') + ^ ^Gq{Q') 
Gm{Q^) = -G2{Q'') , 



(2.42) 



Gq{Q') = Gi(Q2) + G^iQ') + (1 + -^)G3{Q') ■ 

Similarly, one can construct three tensors that project 
out these form factors from the transition amplitude in 
an obvious manner, so that Ge{Q^) = 0^i^^i^A^apiq,q')- 
The amplitude A^af3iq,q') receives contributions from 
many different sources. The electromagnetic proper- 
ties of the vector meson arise from both its quark sub- 
structure as well as from its mixings with other hadron 
states, such as two-pion intermediate states. The quark- 
antiquark contribution to the vector- meson EM form fac- 
tors was explored in Ref. |^] using the same dynamical 
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k+=k+q/2. 
,7'' 




(A, 



~ ----> --" 
k_=-k+q/2 

FIG. 2. Feynman diagram depicting the impulse approxi- 
mate two-pion loop contribution to the p-mes on e lectromag- 
netic transition amplitude A^ai3{p,p') of Eq. (2.43). 



model described in Sec. The EM form factors that 
arise from the quark- antiquark substructure of the vector 
meson are denoted herein as G'f^{Q^). This present study 
extends the work of Ref. by exploring the extent to 
which TT-meson loops contribute to the p-meson EM form 
factors, denoted G^'^{Q'^), and p-meson charge radius. 

In impulse approximation, the two-pion loop contri- 
bution to the £-meson EM transition amplitude is illus- 
trated in Fig. and is given by 



d^k 



whe re k ± = ^q ± k, and the A^*^ are as defined in 
Eq. ( |2.3| ) . Isospin invariance entails that the photon cou- 
ples only to the charged tt meson in the above 7r-loop 
integration. 

The TT-meson EM transition amplitude, given in terms 
of the usual pion EM form factor, is 

r^{k+ + ^Q;Q) = 2{k+ + ^Q)^ FAQ^), (2.44) 

where k++^Q is the average momentum of the incoming 
and outgoing pions, and Q — ~q — q' is the momentum 
transfered to the final vector meso n by the photon 



Substituting the relation in Eq. ( p.44| ) into Eq. ( |2.43[ ) 
one obtains 



X A{k+ + Q)kJ^PP{k_,k+)A{k_) 

X (fc+iQW^^^(fc+ + Q,fc-). (2.45) 

As in the case of the vector meson se lf-ene rgy Ii^^{q), 
the TT-meson loop integration in Eq. ( 2.45| ) necessarily 
samples momenta for which both pions are on-mass-shell 



when q 



ip. It follows that yO- meson EM form 



factors are, in general, complex functions of Q^. 

The imaginary parts of these form factors are associ- 
ated with interference between the photo- stimulated de- 
cay of the p meson into two pions and the usual decay 



of the p meson into two pions. The real parts of the EM 
form factors correspond to the situation where the inter- 
mediate mesons remain off-mass-shell. Hence, they are 
not observed as an asymptotic state and the p meson re- 
mains intact. The real parts of the EM form factors are 
associated with the charge distribution of the p meson. 
Thus, only the real parts of the p-meson electromagnetic 
form factors are of interest to the present study. 

As in the previous calculation of the vector-meson 
self energy 11^^ (q^), a Feynman parametrization is intro- 
duced, thereby combinin g the three denominators of the 
pion propagators in Eq. ( 2.45| ) into a single denominator. 



A(fc+)A(fc+-g)A(fc_) = y daida 



aik-q — {ai + a2){k-Q + jq-Q + \Q^) 



Si 

4 



(2.46) 



The EM form factors of the vector meson are Lorentz 
invariant and can be calculated in any frame. Herein, 
the calculation is performed in the Breit frame for which 
= (0,0, Q,0) and q^ = (0, 0, -iQ, z^iQ2 + ^2)^ 

k ■ Q = xVl^^/Wl^, k ■ q = -\k ■ Q + 



xA(fc+ + Q)A;^-+(fc+ + Q,A:_)A(fc_)A^+-(fc_,A:+), (2.43) 



iz k"^ {jQ'^ + m'j,) , and the elastic condition 2Q ■ q — 
—Q^ is satisfied. Then using the projection operators 
^ua0 discussed above and translating the fc^-integration 



domain by 

k k = k + ^aiq + ^(ai + a2)Q, 



(2.47) 



one finds that the two-pion contribution to the p-meson 
EM form factor is given by 



2i 



(2tt)3 



+1 



k^dk^ \/l--z^dzdxdaida2 
J-i 



[fc2 



k^ 



(2.48) 



where 

k^ - 



-l{l-al)q^-ml + l{ai+a2)ia2 + l)Q^, (2.49) 
x^^.(Q^).9^..r^^(^-,fc+)r^^(fc+ + g,fc-) (2.50) 



)9p7rn 

k' — k-\-^aiq-\-^{ai-\-a2)Q, 
k" = k' + \Q 

fc± = ± [fc + \aiq + i(ai + a2)Q\ 



\q. 



(2.51) 
(2.52) 
(2.53) 



Performing the k integration by parts twice, Eq. (2.4S) 
yields 



Gr(Q') 



(2tt)3 



+1 



\Jl — z'^dz dx dai da 2 



d 



T2"l/c2=0 



(2.54) 
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The real part of this amplitude is of interest, which is ob- 
tained from the principal parts of the integrations. As in 
the calculation of the vector-meson self energy Il^i^{q) 
described in previous sections, calculation of the real 
part of such amplitude requires a knowledge of the off- 
mass-shell transition form factor f^^^ {pi,P2) appear- 
ing i n This is the same form factor that appears in 
Eq. ( 2.25|) for the calculation of the vector- meson self en- 



ergy 11 (g ). The determination of this form factor is 



carried out in Sec. IV, and the results of numerical cal- 
culations of the p-meson EM form factors will be post- 
poned until Sec. M. However, it is important to note that 
the fact that f {pi,P2) appears in the calculations 
of both Il^^{q) and Ge{Q^) provides a self-consistency 
check of the dynamical model for f^^^ {pi,P2)- For a 
dynamical model to be acceptable, it must lead to an off- 
mass-shell transition fo rm fa ctor f ^^^ {pi,P2) that when 



substituted into Eqs. (2.24) and (2.48) gives reasonable 
results for both the p-uj mass splitting and the p-meson 
electromagnetic form factor. This issue is discussed fur- 
ther in Sec. ^ 

Since the p-meson EM form factor receives contribu- 
tions from its quark substructure as well as from meson 
loops, the p meson electric-monopole form factor is writ- 
ten as the sum: 



Ge{Q' 



Gf(Q2) 



GTVTT 
E 



G|^(0) 



GTTTT 
E 



(0) 



(2.55) 



Here, the factors in the denominator ensure the unit 
charge of the p meson Ge{Q) — 1- 

In Ref . 1^ , it was noted that the quark-photon vertex 
obeys the Ward identity which results from the U{1)- 
gauge invariance of QED. When quark-photon vertex sat- 
isfies this identity and the Bethe-Salpeter amplitudes for 
the p meson are normalized in the canonical manner, the 
charge conservation condition [G^''(0) = 1] for the p me- 
son is guaranteed. A similar Ward identity constrains 
the behavior of the pion-photon vertex, 

limr^(p,p-Q)--^A-i(p). (2.56) 
Q^O dp^ 

This identity is satisfied by the form introduced in 
Eq. (2.44). It is straight-forward to show that this iden- 
tity provides a relation between the magnitude of the 
form factor G'^{Q) and the derivative of the self energy 
H^^(— m^, TO^, m^r). This is done by examining the am- 
plitude k^apiq.q') of Eq. ( p. 45 ) in the limit that the 
photon momentum Q — > 0. Then, with the identity in 



Eq. (2.56), one can show that 



^im oJf;'^A^„0(g,Q-q) 



where SZ.^.,^ is that part of the contribution to the nor- 
malization of the p meson that arises from the two-pion 
intermediate state: 



^H [q ,m^,m7r) 



(2.57) 



It follows immediately that G|/^(0) = (J^tttt, which is th e 
justification of the choice of normalizations in Eq. (2.55). 

From Eq. (2.55), one may calculate the p meson EM 
charge radius in the usual manner. 



(r^ 



-6 —,Ge{Q^) 

-6 / d 
l + G|-(0) 



(2.58) 

Gi^(Q') + ^G^^(Q'))-(2.59) 



The charge radius that resul ts fro m a numerical calcula- 
tion of Eq. ( ^.59 ) using Eq. ( ^.54 ) is given in Sec. 0. 



III. APPROXIMATE REDUCTION TO 
TIME-ORDERED CONTRIBUTIONS 

The approach to evaluating the meson-loop corrections 
to the p-uj mass splitting used here is to evaluate the 
covariant Feynman diagram of Fig. which includes 
all possible time orderings of the vertices and interme- 
diate meson propagators. This is in contrast with the 
work of Ref. 0] , which evaluates the mass splitting using 
time-ordered perturbation theory (TOFT), keeping only 
that part of the Feynman diagram corresponding to two 
mesons propagating forward in time (the forward-forward 
time ordering). In TOPT there is a separate contribu- 
tion from the two-meson Z graph, where the vertices are 
ordered in such a way that the initial process (for, say, 
the p self energy due to two-pion loops) is the creation of 
pTTTT from the vacuum, which then propagate forward in 
time with the original p meson until the original p meson 
annihilates with the two tt mesons. This is the backward- 
backward time ordering. A calculation of the forward 
time ordering requires knowledge of the ptttt vertex at 
and near the vicinity of the real decay kinematics, where 
it is reasonably well known. However, a calculation of the 
backward time ordering requires knowledge of the vertex 
for the vacuum to ptttt process, which is not well known. 
Fortunately, the vertex for ptttt creation and annihilation 
is believed to be strongly suppressed, and the backward- 
backward diagram is further suppressed by the energy 
denominators of the ppmr intermediate state propaga- 
tion. For these reasons, the Z-diagram contribution is 
excluded from the calculation described in Ref. 

With the belief that the multiple-hadron creation 
and annihilation vertices should suppress the backward- 
backward time-ordered diagrams, a simple calculation 
is undertaken to assess the relative importance of this 
diagram to the forward-forward time ordering in the 
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present covariant framework. One might argue that pre- 
dictions obtained within this covariant framework would 
be neither robust nor rehable if the forward-forward and 
backward-backward contributions were both large, but of 
opposite sign and canceled each other. Such subtle can- 
cellations would be very sensitive to details of the model 
which may not be constrained accurately enough to be 
trusted. If there were a cancellation between the forward- 
forward and backward-backward diagrams in the present 
framework, it could explain the discrepancy between the 
small vector-meson mass shift obtained herein, and the 
larger mass shifts obtained in the study of Ref. , which 
keeps only the forward- forward contributions. 

To determine the relative importance of the various 
time-ordered cont ributions, one starts with the expres- 
sion in Eq. (2.24), and then rewrites the p-meson self 
energy as a sum of four contributions arising from the 
four possible time orderings of the two pion propagators. 
One of these terms is identified with the contribution to 
the self energy due to the propagation of two forward- 
propagating (positive-energy) pions. The analytic struc- 
ture and physical interpretation of this term, which is 
analogous to that calculated in Ref. and the three 
other terms is discussed in detail. In Sec. VC, a direct 



numerical evaluation of these terms will show that for 
time-like values of the vector-meson four momentum q, 
the dominant contribution to the self energy 
comes from the term identified with the propagation of 
two forward-propagating pions. Therefore, although the 
inclusion of all four time orderings is necessary to main- 
tain the Lorentz covariance of the present framework, 
the additional time-orderings present in covariant expres- 
sions like Eq. (2.24), are noi responsible for the significant 
discrepancy between the mass shifts obtained herein (to 
be presented in Sec. VC) and those obtained by Ref. 
Rather, in the present quantum field theoretic frame- 
work, one observes that the term in vector-meson self 
energy that arises from two forward-propagating pions 
is, in fact, the dominant contribution. 

Before proceeding further, it is necessary to clarify 
some issues concerning the use of quantum field theoretic 
models formulated in Euclidean space. In the present 
study, observables are obtained by analytically continu- 
ing the arguments of the Schwinger functions (Euclidean 
Green functions) into Minkowski space. In principal, 
this is done by first integrating over the internal mo- 
menta so that the Schwinger function depends only on 
Lorentz-scalar products of the external momenta. The 
Schwinger function can then be analytically continued in 
these scalar products into Minkowski space. For exam- 
ple, once the integration over the internal momenta k is 
carried out in Eq. (2.24), the resulting Schwinger func- 
tion Il^^{q^) depends only on the square of the vector- 
meson four momentum q^. The Schwinger function is 
then analytically continued into Minkowski space by let- 



ting q^ In the present study, the use of model 

forms for the elementary Schwinger functions (e.g., quark 
propagators, Bethe-Salpeter amplitudes and meson tran- 
sition vertices) allows the analytic continuation of exter- 
nal momenta to be performed in a straight-forward man- 
ner by letting external momenta become complex before 
the internal momenta are integrated over. 

An important ramification of such a framework is that 
a direct comparison between the elementary Schwinger 
functions employed herein and those employed in other 
models formulated in Minkowski space is impossible. One 
reason for this is that the model form of the quark prop- 
agators Sf{k), given by Eqs. (4.5), are parametrized in 
terms of entire functions of k. This parametrization 
provides the model with quark confinement by ensur- 
ing that the quark propagator Sf{k) is free of singular- 
ities in the finite, complex-fc plane. Hence, the quark 
propagator has no Lehman representation and no quark- 
production thresholds can arise from the calculation of 
observables. Of course, the use of such quark propa- 
gators also forbids the possibility of performing a Wick 
rotation in the loop-integration variable k. Therefore, 
the elementary Schwinger functions employed herein are 
permanently embedded in Euclidean space. Only final 
quantities, such as 11^^ (q^), can be analytically contin- 
ued into Minkowski space. 

Another feature of the framework that makes a Wick 
rotation impossible to perform is the rich analytic struc- 
ture of the meson transition vertices, e.g., A'^^'^{p,p'). A 
Wick rotation of these transition vertices, although possi- 
ble in principal, is made difficult by the presence of singu- 
larities which arise from the many-particle nature of the 
field theory. In general, these singularities may cross the 
integration path requiring that exceptional care be taken 
when performing a Wick rotation. In practice, these 
transition vertices are calculated and parametrized only 
on the small domain of the complex-momentum plane 
for which they are needed in pion-loop integrations; their 
behavior outside this domain is unknown. Such difficul- 
ties are not present in a quantum mechanical formulation 
that truncates the Fock space in such a way as to include 
only states with a small number of particles. The an- 
alytic properties of transition vertices obtained in such 
models are more easily analyzed. 

Nevertheless, it is possible to make a simple connection 
between this field theoretic approach and TOPT. The 
starting point is the calculation of th e p-n ieson self energy 
due to the two-pion loop from Eq. ( 2.24 ). As in Ref. Q, 
the calculation is carried out in the rest frame of the 
p meson. In this frame, = {0,imp) and the 7r-meson 
propagators, which depend on the momenta k± ~ fc± ig. 



* Additional details concerning the formulation of quantum 
field theory in Euclidean space and the analytic properties of 
Schwinger functions may be found in Ref. |13|. 
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can be written in the following form 
A(fc±,m^) = 



(3.1) 



This nomenclature specifies the four possible time order- 
ings of the two pions in the pion loops. It is justified from 
an analysis of the motion of the poles in the complex-fc4 
plane as the three-momentum k is varied. 



where Lj{k^) — yk'^ + m^. The denominator of Eq. (3.1) 
is the difference of two squares. It can be factorized and 
rewritten as a difference of two poles in the complex-fc4 
plane, 



Im(k4 ) 



2m, 



A(A;±,m^) = 



1 



2w(/c2 



1 



iki ^ \mp + uj{k^) — it 



ik4 =p ^iTip — uj(k'^ 



(3.2) 



In a field th eory based in Minkowski space, the two 
terms in Eq. (3.2) would be identified with forward- and 
backward-propagating pions. A similar identification is 
possible for the two-pion syste m in our approach by sub- 
stituting Eq. (3_^) into Eq. ( 2.24 ). One then obtains 
an expression for the vector meson self energy in which 
four poles in the complex- A:4 plane are explicit. Writ- 
ing the self energy as a sum of four terms: Il^^{q'^) — 
n[++l(g2) + n[— I(g2) + n[+-l(g2) + n[-+l(92), where, for 
a.b = +, — . one has 



X 



X 



X 



X 



.Re(k,) 



Imik4}^ 
2m^ - 

nil - 

x\ 



Y 



,Re(k, 



1 — -m. 



1 — -2m„ 



X 

A 

FIG. 3. The positions in the compIex-fc4 plane of the poles 



in Eq. (3.3). Dashed lines with arrows indicate the motion of 



^- P' 3 J (27r 



d-^k T^^{k,q) 



4w2(fc2) 



and 



A[++l(fc) 



A^—\k) 



iki — ^rup + cj(fc2) — ie I 

V 

yik4 + ^nip — uj{k'^) + it 

• . ) 

iki — \'>Tip — Lu{k'^) + it J 



ik4 + ^rrip + Lo{k'^) — it 



these poles as the relative three-momentum \k\ is decreased 
from infinity to zero. The positions of the poles for |fe| = are 
denoted by "X." The diagram on the left depicts the fictional 
A ["''1 ( fc ) , (3.3) situation for which nip « < 2m,r and the p — > tttt decay 

is forbidden by energy conservation. The diagram on the right 
depicts the situation for which mp > 2m,r and the p — > tttt 
decay is allowed. For a certain values of \k\, two of the pion 
\ poles may pinch the fc4-contour integration as e ^ 0. This 

signals the threshold of the p — > tttt decay and results in a cut 
in the self energy of the p meson. These two poles correspond 
to the forward- forward time ordering in Eq. (3.4). 



(3.4) 



(3.5) 



A^+-\k) 




A^-+\k) 



ik4 + \rap + uj{k'^) — it 

- ' , ) 

ik4 — ^irip — ULi[k'^) + it I 

( 

\ik4 + \mp — w(fc2) -I- 



(3.6) 



(3.7) 



These are identified as contributions to the p-meson 
self energy due to forward- forward (++), backward- 
backward ( ), forward-backward (H — ) and backward- 
forward (— +) propagation of the two pions, respectively. 



These integrands in Eq. (3.3) exhibit poles for particu- 
lar values of the fourth component of the loop-integration 
momentum ^4. The positions of these poles depend on 
the magnitude of the three-momentum k as well as the 
mass of the vector meson mp. In Fig. ^, the positions of 

these poles, along with their motions as |fc| is varied, are 
depicted for two possible scenarios. The left diagram cor- 
responds to a scenario (not realized in nature) for which 
the p-meson mass nip « rriT^ /2 < 2to^. In this situation, 
the center-of-momentum energy of the p meson is below 
the two-pion decay threshold; hence, the p meson is sta- 
ble with respect to the strong decay into two pions. The 
diagram on the right depicts the situation (realized in na- 
ture) for which the mass of the p meson nip > 2mT^ and 
the strong, p — > tttt decay is allowed. In both diagrams, 
the motions of these poles as is varied are denoted by 
dashed lines. Arrows depict the trajectory of the poles 
as |fc| is decreased from infinity to zero. The positions of 
the poles for |fc| = are denoted by "X." 

Examination of the left diagram in Fig. depicts the 
fictional situation where the p-meson mass is below the 
p —t nn threshold energy. Hence, the integrations over 
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ki in Eq. ( ]3.3| ) never encounter poles for any value of 
Thus, in this scenario, the p-meson self energy 11^^ {q'^) 
is purely real and the p meson is stable with respect to 
strong interactions. In the right diagram of Fig. ^ for 
which TUp > 2m7r, one observes a very different situation. 
Here, as the relative momentum |fc| is decreased from 
— oo, two of the poles pass each other, thereby causing 
the contour integration path for to become "pinched" 
in the limit that e ^ 0. Such a pinching of the k^- 
integration contour results in a branch cut in 11^^ (g^) 
and hence a non-zero imaginary part of the vector-meson 
self energy. This signals the opening of the two-pion de- 
cay channel and provides the p meson with a finite two- 
pion decay width. The value of |fc| for which this pinch 
first develops is obtained from the energy-conservation 
condition: uj{k'^) = mp/2. The fact that this pinch sin- 
gularity in n-^^(g^) can only arise in the calculation of 
n++(g^) is the justification for this term being identi- 
fied with the propagation of two forward-propagating pi- 
ous. Clearly, only positive-energy pious may appear in 
asymptotically free states and in so doing contribute to 
the imaginary part of the vector-meson self energy. 

The term 11 {q'^) is identified as arising from two- 
backward propagating pious by noting that had one 
started with the convention that energies were negative 
rather than positive quantities, the replacement of the p- 
meson center-of-momentum energy nip with —mp, would 
result in the pinch singularity appearing in 11 (q^) 
rather than in n++((7^). That is, within a framework 
that employed such a convention, the asymptotically free 
states would necessarily contain negative-energy pious. 
A corollary that follows from the above argument is that 
in the limit that q^ 0, no energy flows into the two- 
pion-loop diagram and there is no difference between the 
forward-forward and backward-backward time orderings; 
that is, for a massless vector meson these time orderings 
contribute equally to its self energy. (This is indeed what 
is observed numerically, as discussed in Sec. VC.) 

It is important to note that the contribution to the 
vector-meson self energy due to negative energy pious 
does indeed appear in quantum mechanics, even though 
there are no negative energy states in the Fock space. 
By considering the time-ordered Feynman diagram asso- 
ciated with n~~((7^), written so that all particles move 
in the positive-time direction, one observes that the four- 
particle Fock state Ippmr) has a non-zero overlap with 
the intermediate particles depicted in the diagram. This 
state contributes only to the real part of the vector-meson 
self energy because energy conservation forbids the decay 

p —^ ppTTTT. 

In a quantum mechanical framework, such as that em- 
ployed in Ref. Cauchy's Theorem entails that the 
mixed, forward-backward (-1 — ) and backward-forward 
( — h) contributions, are identically zero. However, in 
the present framework the meson-transition vertex in 
Eq. (4.9) and quark propagator in Eq. (4.5) have essential 



contour the integration at infinity and use Cauchy's the- 
orem to show that these terms are zero. Consequently, 
there is no such constraint on the value of these terms in 
the present approach, and they will generally be different 
from zero. A discussion of the physical interpretation of 
these terms and the role they play in determining the 
unitarity of the theory is left for future investigations. 
It is sufficient for the present application to note that 
in Sec. VC, a direct numerical evaluation of these terms 



singularities at infinity. Therefore, one cannot close the 



shows them to provide a negligible contribution to the 
vector-meson self energy. Hence, they have no impact on 
the results presented herein. 



IV. MESON TRANSITION AMPLITUDES 

In the previous sections, expressions for the vector- 
meson self energy and EM form factors were given in 
terms of meson propagators and meson transition ampli- 
tudes. In the following section, these transition ampli- 
tudes are given in terms of nonperturbative quark-loop 
integrations that depend on model Bethe-Salpeter ampli- 
tudes and nonperturbatively-dressed quark propagators 
developed in studies of hadron observables based on the 
Schwinger-Dyson equations of QCD. A recent review of 
such studies is provided by Ref. fl^ . 

The imaginary part of the vector-meson self energy is 
determined by meson transition amplitudes for which all 
external meson momenta are on-mass-shell. It follows 
that an experimental determination of decay widths for 
the vector meson provides only the on-mass-shell values 
of the meson transition amplitudes. In the following, the 
values of the transition amplitudes when all mesons are 
on- mass-shell are referred to as coupling constants, e.g., 

A determination of the real parts of the vector meson 
self energy requir es ca lculati ng th e principal part of the 
integrals in Eqs. (2.24) and ( 2.28| ), which sample the k^- 
integration domain for which the intermediate mesons are 
off-mass-shell. Thus, a calculation of the real parts of the 
self energy and EM form factors requires knowledge of the 
meson transition form factors for off-mass-shell values of 
the intermediate meson momenta. As the real part of the 
self energy is not directly observable, experiment cannot 
determine the behavior of the off-mass-shell-meson tran- 
sition form factors. Nonetheless, experimental observa- 
tions can provide some guidance in helping to constrain 
the real parts of self energies and EM form factors. This 
can be done by looking for observables which are propor- 
tional to th e diff erence of tw o self energies. For example, 
from Eqs. ( 2.30D and (2.31), it is apparent that the p-Lo 
mass splitting provides a measure of the difference of the 
real parts of the p and lo self energies. In Sec. IV , the cal- 
culation of off-mass-shell meson transition form factors is 
discussed. 

The VVP and VPP vector-meson transition ampli- 
tudes are calculated within a generalized impulse ap- 
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proximation in which the mesons couple to each other 
by means of a nonperturbatively-dressed quark loop. In 
this generalized impulse approximation, the quark propa- 
gators and meson Bethe-Salpeter amplitudes are nonper- 
turbatively dressed, but explicit three-body interactions 
are neglected. 

In the following, quark- loop integrations are performed 
in Euclidean space (where > 0), and final results are 
analytically continued to Minkowski space by letting the 



?/. The Dirac ma- 

t - 



external (timelike) momenta q —m 
trices 7^ employed herein are Hermitian, 7^ — 7/j, and 
satisfy the anti-commutation relation {7^^,71^} — ^Sfj^i^. 

The amplitude describing the coupling of a vector me- 
son of momentum q = —pi — p2 and two pseudoscalar 
mesons with momentum pi and p2 is given in impulse 
approximation by^, 

j^S{k++)X,Vf,{k+^p2;q) 

xS(fc_+)A,rp(fc-ig;pi)S(fc__)Aferp(fc;p2), (4.1) 

where kafs — k + ^q + ^P2, K are linear combinations of 
the Gell-Mann SU(3)-flavor matrices associated with the 
mesons involved, the trace is over color, flavor and Dirac 
indices, and S(A:) = diag(S'tj(fc), Sd{k), Ss{k)) is the quark 
propagator for u, d and s quarks. It is a 3x3 matrix in 
the quark-flavor space. In this study, the effects of SU(3)- 
flavor breaking are neglected, so the quark propagator is 
diagonal in flavor indices and Su{k) — Sd{k) — Ss{k). 
Here V^(/e, k') is the vector- meson Bethe-Salpeter ampli- 
tude for the p or a; mesons, and rp(fc;pi) is the Bethe- 
Salpeter amplitude for the P = n or K mesons. The 
flavor structure of the me son Bethe-Salpeter amplitudes 
is given explicitly in Eq. (O) in terms of the 3x3 Gell- 
Mann matrices in the quark-flavor space. The meson 
transition amplitude A^'^(pi,p2) is depicted in Fig. |l] as 
a triangle diagram. 

The transition form factor defined by Eq. ( ^.3| ) is ob- 
tained from Eq. (4.1) by choosing values of the flavor 
indices i, j and k appropriate for the particular process 
V PP under consideration, and then contracting the 
amplitude with an appropriate vector. In the case of 



one obtains 



[Pl + P2 



^{puP2), (4.2) 



where the flavor indices refer to the following linear com- 
binations of Gell-Mann matrices A± = -^(Ai ± 1X2)1 
which are associated with an isovector meson of good 
charge. The resulting form factor f^^^ {pi,P2) is the 
essential dynamical element necessary to calculate the 



^The sense of all momenta is such that positive momenta 
flow into vertices. 



real part of the two-pseudoscalar intermediate state con- 
tribution to the vector- meson self energy Il^^{q^) from 
Eq. (|^. 

Transition amplitudes for processes involving strange 
mesons, such as p — > KK, are obtained by assuming 
SU(3)-flavor symmetry is exact, so that Ss{k) — Su{k) 
and rj<-(fc;pi) = Tjrik^pi). From considerations of 
SU(3)-flavor symmetry, one can show that the left-hand- 
side of Eq. ( |4.2| ) gains an additional factor of l/\/2 for 
the process p — > KK. In p artic ular, for p+ K^K^ , 
the analogous result for Eq. (4.2) would be 



-^I (Pi,P2) - — ZTT^M (Pi^P2)- (4.3) 



V 



Here, the flavor indices denote the linear combinations 
of Gell-Mann matrices, A^- = -^[X^ — iX^) and X^o = 

-^{XQ+iXi). The additional factor of l/\/2 that appears 

in Eq. (^^) relative to that in Eq. ( |4.2| ), is the reason 
for the factor of 1/2 that appears as the coefficient of 
n^^((7^mK,mA') in Eqs. (|]30|) and ( |2.3lD . 

Similarly, the coupling of two vector mesons of mo- 
menta — (pi +P2) and Pl to a pseudoscalar meson of mo- 
mentum p2 is given in the impulse approximation by 



A;i^(Pi. 



P2) 



/d'^k 
^^;^S(fc++)AjV^(fc+ip2;'?) 



xS(fc_+)A,K(fc-ig;pi)S(fc__)Aferp(fc;p2). (4.4) 



The eleinents appearing in Eq. (4.4) are the same as 
in Eq. (4T). No new elements are necessary to ob- 
tain the VVP transition amplitude^ ^The form factor 

/^^■^(Pi,P2) is obtained from Eq. ( ^!^ using Eq. ( |2^ ) 



and is the essential dynamical element necessary to 
calculate the real part of the intermediate vector- 
pseudoscalar contribution to the v ector-meson self energy 
H^^(g2,my,mp) from Eq. ( |2^ ). 

Here too, the assumption of exact-SU(3)-flavor invari- 
ance provides a simple relation between the processes 
p — > Lun, p — > prj, p — !■ K*K, lo — s- p7r, lu — s- Lurj, and 
uj — !■ K*K. The flavor factors that relate these processes 
lead to the coefficients appearing i n fron t of t he ve ctor- 
pseudoscalar contributions in Eqs. ( ^.30| ) and ( 2.31 ) 



Evaluation of the quark- loop integrations in Eqs. (4.1) 
and (4.4) requires knowledge of the m-, d- and s- 
quark propagators as well as the vector-meson and 
pseudoscalar-meson Bethe-Salpeter amplitudes. These 
are taken from phenomenological studies of the electro- 
magnetic form factors and strong and weak decays of 
pseudoscalar and vector mesons [^[H . 

The form of the model, dressed-quark propagators 
are based on several numerical studies of the quark 
Schwinger-Dyson equation using a realistic model-gluon 
propagator |10 1^. In a general covariant gauge, the 
propagator for a quark of flavor / is written as Sf{k) = 



—17 • kay{k ) + <Tg{k ). It is well described by the fol- 
lowing parametrization: 
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TABLE I. Confined-quark propagator parameters for u, d, 
and s quarks from Ref. An entry of "same" for a meson 
indicates that the value of the parameter is the same as for 
the quark flavor with which it is grouped. 



//meson 


"o 


t{ 




bi 


Cf 


nif [MeV] 


u, d 


0.131 


2.900 


0.603 


0.185 





5.1 


TV 


same 


same 


same 


same 


0.121 





p, u 


0.044 


0.580 


same 


0.462 








s 


0.105 


2.900 


0.540 


0.185 





130.0 


K 


0.322 


same 


same 


same 








K* 


0.107 


0.870 


same 


0.092 









+2fnf£,[2{x + mj)] + Cje" 
2{x + fh)) - 1 + e-^(="+™/) 



2\2 



2(.T + TOp 



(4.5) 



where £^[x\ — (1 
rhf 



-J 



4 



= ruf/X, A = 0.566 GeV, A = 10 

and where the b{ , Cf and m/ are parameters. Since 
£_[x] is an entire function, this dressed-quark propagator 
has no Lehmann representation. The lack of a Lehmann 
representation for Sf{k) is sufficient to ensure that the 
quarks are confined, since no quark-production thresh- 
olds can appear in the calculation of observables. The 
quark propagator Sf{k) also reduces to a bare-fermion 
propagator when its momentum is large and spacelike, in 
accordance with the asymptotic behavior expected from 
perturbative QCD (up to logarithmic corrections). The 
parameters for the d- and s-quark propagators were 
determined in Ref. Q by performing a fit to ^ range of 
TT- and X-meson observables. The parameters for the s- 
quark propagator employed herein were re-fit in Ref. , 
also from a fit to iiT-meson observables, but with the 
additional constraint that the resulting s-quark propaga- 
tor more closely resemble that obtained by the numerical 
study of Ref. ||ll| . The resulting parameters are given in 
Table |. 

In the chiral limit m/ — > 0, Goldstone's theorem allows 
one to obtain the pseudoscalar-meson Bethe-Salpeter am- 
plitude Tp{k;P) directly from the dressed-quark propa- 
gator Sf{k) [|ll|. For the calculation of infrared transi- 
tion amplitudes, like those studied herein, it is sufficient 
to keep only the part of the pseudoscalar Bethe-Salpeter 
amplitude that is proportional to 75. In this case, the 
amplitude is given by 



rp(fc;P)=75 



Bpjk'^ 
Np 



(4.6) 



where Bp{k^) has the same form as Bf{k^), one of the 
Lorentz invariants appearing in the quark inverse prop- 
agator Sj^{k) = Z7 • kAfik"^) + Bflk"^), but with the 



parameters given in Table | for P = n or K . The tt- and 
if-meson Bcthc-Salpeter amplitudes are obtained from 
Eqs. (4.5) and (4.6) using the parameters from Table |. 

In Eq. (16), Np is the on- mass-shell normalization 
for the P = TT or K meson, and is determined by the 
normalization condition: 



1 



mp 



tr 



D 



d^fc 



— S{k+y)Tp(k; p)S{k- y)fp{k- p) 
+S{k+y)Tpik;p)^S{k~^p)Tpik;p), 



(4.7) 



where tr^; denotes a trace over Dirac indices and Nc — 3. 
In Ref. [ pT| it was shown that by maintaining only the 
Dirac amplitude proportional to 75 in the pion Bethe- 
Salpeter amplitude, one cannot reproduce the ultravio- 
let behavior of the pion electromagnetic form factor. To 
do so requires the inclusion of the axial- vector 757 • p 
components into the 7r-meson Bethe-Salpeter amplitude. 
Nonetheless, for the calculation of observables which in- 
volve small momentum transfers, neglecting Dirac mo- 
ments other than 75 and nor mal izing the Bethe-Salpeter 
amplitude according to Eq. (4.7) is usually sufhcient. 

The vector meson Bethe-Salpeter amplitudes employed 
herein are taken from Ref. Q : 



V^{k;p) = 7m + 



Pfil-P\ Byjk^) 
Nv 



(4.8) 



where p is the momentum of the vector meson with 
mass my, k is the relative momenta of the quark and 
antiquark, Ny is the Bethe-Salpeter normaliz atio n, ob- 
tained from an obvious generalization of Eq. (4.7), and 
Bv{k^) has the same form as the Lorcntz-invariant func- 
tion Bf{k^) found from the quark propagator of Eq. (4.5) 
with parameters given in Table || with V — p, oj or K* . 

Again, all but the leading Dirac moment proportional 
to 7^ in the vector meson Bethe-Salpeter amplitude have 
been neglected. A numerical study of the Bethe- 
Salpeter equation that employed a separable model for 
the quark-antiquark scattering kernel found that the 
magnitude of this moment is about ten times larger than 
the next largest. However, some observables may be more 
sensitive to sub-leading Dirac moments in the Bethe- 
Salpeter amplitude. Indeed, this was the case observed 
in Ref. for the asymptotic-g^ behavior of the tt meson 
electromagnetic form factor. However, the simple form 
provided by Eq. ( |4.8D is sufficient to provide a reasonable 
model for calculating the off-mass-shell transition form 
factors /^^^(pi,P2) and /^^^(pi,p2). 

With these elements having already been deter- 
mined by previous studies, the on- and off-mass-shell 
values of the transition amplitudes A^^^ {pi,p2) and 
A^y^{pi,P2) can be calculated with no adjustable pa- 
rameters. The value of a transition amplitude for on- 
mass-shell meson momenta is referred to as a "coupling 
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constant" . If the invariant mass of the final state ^/s is 
less than the mass of the initial state, the decay is allowed 
to proceed and the coupling constant is experimentally 
observable. If energy conservation prevents all external 
mesons from being simultaneously on-mass-shell, one can 
still define the "coupling constant" as the value of the 
transition amplitude at this point. However, in this case, 
there is no decay from which it can be measured experi- 
mentally. 

In Ref. H], the coupling constants (?p7r7r and gx'KTr 
were obt ained using the model quark propagators from 
Eqs. (|4.5| ) and Bcthe-Salpeter amplitudes from Eqs. ( |4.6| ) 
and ( |4.8| ) with the parameters from Table |. The values 
for QpT^TT and qk'K-k obtained therein were 8.52 and 9.66, 
respectively. The experimentally observed values for 
these coupling constants are 6.03 ± 0.02 and 6.40 ± 0.04, 
respectively. There are several possible reasons that this 
model gives results for the decay widths of a vector me- 
son into two pseudoscalars that are larger than obtained 
by experiment. The first is that for this process in par- 
ticular, higher Dirac moments contribute significantly to 
the magnitude of the coupling constant ||T^ . The second 
is the possibility that final-state interactions may also be 
important. For these reasons, rather than employ the 
value of gpT^TT obtained in Ref. § , the experimentally de- 
termined value of gpTTTT = 6.03 is used. The meson transi- 
tion form factors f^^^ {pi,P2) and f^^^ {pi,P2) are the 
important dynamical element necessary for the present 
calculations. The scales of these are rather insensitive to 
details of the model Bethe-Salpeter amplitudes, so that 
the model of Ref. H] is quite adequate for their determi- 
nation. 

The coupling constant gp^jTr = 10.81 is obtained from 
the pLun transition amplitude in Eq. (4.4), when all three 
mesons are on-mass-shell. Of course, the lack of phase 
space makes the decay to pir impossible, so that 
the vector-meson self energy contribution due to vector- 
pseudoscalar intermediate meson states is real for a phys- 
ical vector meson; i.e., the value of gpunr cannot be de- 
termined experimentally. Throughout this articl e, the 
value 5p^^ = 10.81 obtained from Eqs. ( |2.28D with ( |2.29| ) 
is used . For comparison, the value of this coupling 
constant can been estimated from the anomalous Wess- 
Zumino term in an effective chiral Lagrangian with the 
assumption of a universal p-meson coupling. With the 
definition of the coupling constant from Eq. (2.8), it is 
found to be ^p^^r ~ 11.5 ||l^, which is consistent with the 
results of the present calculation. Similarly, a study of 
the decays of the cj-meson using vector-meson dominance 
finds the value gpuj-a- « 9.7 Ref. 

The loop integrations of Eqs. (2.24) and (pT28) sam- 
ple the meson transition form factors {pi,P2) 
and f^^^ {pi,P2) for off-mass-shell values of the 
intermediate-meson four momenta pi and p2- A calcu- 
lation of thes e for m fac tors from the quark-loop integra- 
tions of Eqs. (4.1) and (4.4) requires a suitable definition 
for meson Bethe-Salpeter amplitudes when their total 
momentum is off-mass-shell. Although Bethe-Salpeter 
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FIG. 4. The me son transition form factor /^^^(pi,P2) 
obtained from Eqs. ( [l.l[ ) (circles) and the numerical fit (solid 
curve) of Eq. ( [l.ll[ ) used in the calculations of vector meson 
self energies and EM form factors. 



amplitudes are uniquely determined from the homoge- 
neous Bethe-Salpeter equation for on-mass-shell values 
of the total momentum, an extrapolation of the ampli- 
tude away from this on-mass-shell value is not uniquely 
determined, and hence, is model dependent. 

Her ein , the Bethe-Salpeter amplitudes as defined in 
Eqs. (4.6) and ( |4.^ ) are used, even for values of the me- 
son four momenta that are off-mass-shell. This is most 
similar to the procedure employed in quark model calcu- 
lations such as Ref. Q , where a complete set of interme- 
diate states is inserted between the vector-meson states 
leading to transition amplitudes for on-shell mesons. The 
off-shell behavior of the intermediate mesons is then de- 
termined solely by the behavior of the meson propaga- 
tors. 

Rather than directly use the numerical values of 
the meson transition form factors as calculated from 



Eqs. (4.1) and (4.4), it is found that the resulting form 
factors can be parametrized in terms of the following 
hmctions: 



exp 



fc2 



9 19 

^vppy 



1 — TOp/6' 



(4.9) 



(4.10) 



where pi ~ + k and p2 ^ — k in Eq. (4.9) and 
Eq. (4.10). The parameters are chosen by fitting these 
forms to the numerical results on the domain fc^ € (0, 
1) GeV^ for a range of external vector- meson four mo- 
menta q"^ G (—1, -1-1) GeV^. The g^-dependence of the 
resulting parameters is well described in terms of the lin- 
ear forms: 
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FIG. 5. The meson transition form factor / obtained 
from E qs. ( 1.4) (circles) and the numerical fit (solid curve) 
of Eq. (4.12) used in the calculations of vector meson self 
energies. 



bvpp{q^) = feypp + bypp 
h^lpp = 0.855 GeV, 

a 

'VPP 



6lpp = 0.0693 GeV-\ 



(4.11) 



and 



bvvp{q^) = bvvp + bvvp q^ 
b°yvp = 1-098 GeV, 
h\^yp = 0.159 GeV-^ 



(4.12) 



The numerical results (denoted by circles) and the fits 
(solid curves) are shown in Figs. ^ and ^ for = —m^ 

for the form factors f^^^ipi,P2) and f^^^{pi,P2), re- 
spectively. The similarity of the two curves suggests that 
both could have just as easily been fitted to either expo- 
nential or dipole forms. The choice of one form over the 
other is arbitrary, and does not affect the results obtained 
herein. 



V. RESULTS AND DISCUSSION 

A. mass dependence 

Fig. 1^ shows the real (solid curves) and imaginary 
(dashed curves) parts of the contributions to the vector 
self energy due to tttt and KK loops as a function of the 
four- momentum squared of the vector meson . These 
curves are obtained by calcul ating Il^^{q ^, m,p , mp) 
with P = TT or is: in Eqs. ( |2.24D and ( ^.25] ) with 
the assumption of exact SU(3)-flavor symmetry; i.e., 
(jfp^^ = gpKK = 6.03 and using the parametrization of 
r''''ipi.P2) given by Eq. (|I|). 



FIG. 6. The real (solid curves) and imaginary (dashed 
curves) parts of the vector-meson self energy due to inter- 
mediate TTTT and KK meson states. The vertical dashed line 
denotes the on-mass-shell point — --m^p for the p meson 
and the vertical solid line denotes q^ = 0. 



For spacelike values of the vector-meson momenta 
> (to the right of the solid vertical line in Fig. |^), 
the self energy 11^^ (g^, mp, top) is purely real and ap- 
proaches zero monotonically from below as q^ +oo. 
The threshold for the decay p ^ tttt is at q^ = — 4toJ 
(just left of the solid vertical line in Fig. ||). Here, the 
imaginary part of the self energy n^^(g^, to^, TO,r) be- 
comes non-zero. This abrupt change of the imaginary 
part of the self energy at the two-pion threshold leads to 
a discontinuity in the second derivative of the real part 
of the self energy. This causes the real part of the self 
energy to rapidly turn over at q^ w —0.3 GeV^ and start 
to approach zero as q^ becomes more negative (time- 
like). The same behavior is observed in the two-JsT-meson 
self energy Il^^{q^,mK,mK), but much deeper in the 
timelike region since the two-kaon decay threshold is at 



q 



1 GeV2. 

The discontinuity in the second derivative of 
Ren^^((7^, TOp, TOp) at q^ = —Anip is a general feature 
of an amplitude near the threshold of a decay into a two- 
meson intermediate state with L = 1 relative angular 
momentum. In general, the dependence of the imaginary 
part of the self energy on the center-of- momentum energy 

^/s near a decay threshold is proportional to -^/s^^^^^^ 
where L is the relative angular momentum of the two 
outgoing decay products. The centrifugal barrier due to 
the relative angular momentum L tends to soften the be- 
havior of the imaginary part of the amplitude at thresh- 
old. In the present case of p — > tttt, the pions have 
relative L = 1 angular momentum and hence the sec- 
ond derivative of Ren^^(g'^, top, mp) is discontinuous at 



-Am^p. 



Had this two-body decay process proceeded 
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0.0 

q (GeV') 

FIG. 7. The real (solid curves) and imaginary (dashed 
curves) parts of the vector-meson self energy due to interme- 
diate ijj-K and K* K meson states. The vertical dashed line 
denotes the on-mass-shell point - 
and the vertical solid line denotes 



—nip for the p meson 
0. 



in the s-channcl with L = 0, the first derivative of the self 
energy would have been observed as being discontinuous 
at threshold = — 4mp. 

In Fig. the real (solid curves) and imaginary (dashed 
curves) parts of the vector-meson self energy due to 
vector-pseudoscalar intermediate state 11^^ (q^, mv,mp) 
is shown as a function of q^. These curves are ob- 
tained by calculating 11^^ (g^, my, mp) with P — t t and 
V^poTP^K£indV^K*m Eqs. ( ^^ and ( |2^ ) 
with the assumption of exact SU(3)-flavor symmetry; i.e., 
gpujTT = gpK'K = 10.81 and using the parametrization of 
/^^^ (pi , P2 ) given by Eq. . 

Although not apparent from Fig. 0, when the 
vector-meson four momentum becomes large and space- 
like, q^ +00, the real part of the self energy 
n^^(g^, TOy , mp) approaches zero from above. The fact 
that the decay threshold for p lutt is observed at 
q^ = —{moj + m^)^ w -0.84 GeV^, so far above the 
on-mass-shell value of q^ — —mj, means that very lit- 
tle structure for the real part of n^^((7^, TOj^, m^) is 
observed for the range of q^ plotted in Fig. ^ The 
lack of structure in the real part of li^^{q^,mK*,TnK) 
for the K* K channel is even more apparent, owing to 
the fact that the threshold for p K*K is found at 
g2 = -{mK* +mK)^ w -1.9 GeV^ (not shown in Fig. 0). 

As discussed above for 11^^ (q^ ,mp,mp), the order in 
which the discontinuity in derivatives of the real part of 
the self energy appears is determined by the relative an- 
gular momentum L of the out-going decay products. It is 
straight-forward to show that the decay of a vector me- 
son with = 1~ into a vector and pseudoscalar meson 
with = 0^ must be in an orbital angular momen- 



tum L = I state. Hence, the discontinuity is expected to 
appear in the second derivative of Ren^^(g^, my, ?Tip) 
at q^ — ~{m^ -|- m^)^. Numerical differentiation of the 
results shown in Fig. |^ confirm the presence of this dis- 
continuity. 

One observes that 11^-^ (g^, my , mp) passes through 
zero at = 0, a feature that is not observed for 
\lPP{q^,mp,mp). This zero is a result of the structure of 
the VVP meson transition amplitude A^j^(pi,p2) given 
in Eq. (^^), which is required by Lorentz covariance. The 
requirement that this amplitude be proportional to 



f-tivaP Pla P2f3 = £p,}^af3 9q Plf3, 



(5.1) 



leads to a vector-meson self energy 11^^ (g^, my, mp) 
that is necessarily transverse to g^ for all values of q^ and 
hence is of the form given in Eq. (|]^). It follows that 
the vector-meson self energy in this channel vanishes at 
g2 = 0. 

Because the real part of n^jf(g^, my, mp) goes 
through zero at g^ = 0, approaches zero in the deep 
spacelike region and is infinitely differentiable below the 
breakup threshold at g^ = —(my + mp)^, it is a small 
and slowly varying function for g^ > — m^. Therefore, 
the WTT*' intermediate state is expected to contribute very 
little to the p meson self energy (from Fig. |^, it is clear 
that the K*K state contributes even less). However, this 
channel is important for the self energy of the w meson, 
owing to the fact that there are three different pir chan- 
nels that contribute to the w-meson self energy; they are 
p+TT^, /0^7r+, and p'^ir^. This leads to the overall factor 
of three in front of U^^{q'^,inp, m^) in Eq. ( ^.3l| ), which 
makes the pir channel as important for the w-meson self 
energy as the mr channel is for the p-meson self energy 
(see Table ^. 

A single coupling constant and meson-transition form 
factors has been used for each of the two types of inter- 
mediate states {PP and VP) considered herein. The 
result is that the only difference between the various 
channels within a particular type of meson loop is the 
masses of the intermediate mesons involved. A com- 
parison of the self energies resulting from the KK and 
K*K channels with tttt and lutt channels, respectively, 
suggests that states with heavier mesons contribute less 
than states of lighter mesons. For example, for all g^ > 
— m^, the value of 11^^ {q^ ,mK ,mK) is less than half 
the value of H^^(g^, m^r, m^r). The same is true when 
11^^ [q^ ,mK* jItik) is compared with H^^(g^, m^^, m^) 
in Fig. |. 

Within the present model, the dependence of the 
vector-meson self energy on the mass of the intermedi- 
ate mesons can be explored. This is done by letting one 
or both of the masses mp in Eqs. (2.24) and (|2.2j) vary 



while keeping g 



-mj,. The solid curve in Fig. M is the 



result of a calculation of the real part of the self energy 
n^^(— m^, m^, mi) in whi ch on e of the pseudoscalar- 
meson propagators in Eq. (2.24) has mass m,r and the 
other propagator has a mass mi which is varied from 
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FIG. 8. The dependence of the real part of the vec- 
tor-meson self energy on the mass mi of the intermedi- 
ate state mesons. The dashed curve shows the dependence 
of the vector-pseudoscalar channel 11 {—mp,mp,mi) when 
the mass mi of the pseudoscalar meson is varied. The 
solid curve shows the dependence of the two-pseudoscalar 
channel n^^(— m^, m^, mi) when the mass mi of one 
of the pseudoscalar mesons is varied. The short-dashed 
curve shows the dependence of the two-pseudoscalar channel 
— mp,mi,mi) when the mass mi of both pseudoscalar 
mesons is varied. 



to 1500 MeV. It is clear from Fig. ^ that the magnitude 
of the self energy depends critically on the mass mi. In 
particular, for mi = rrip, the real part of the self energy 
— m^,m^,mi) has already decreased by a factor of 
two from its value for mi = m,r- The short-dashed curve 
in Fig. H is the self energy n-^^(— m^, mi, mi) obtained 
with the masses of the two pseudoscalar mesons are taken 
to be equal and varied together. Increasing both masses 
of the intermediate mesons at the same time causes the 
self energy to decrease very rapidly. 

It was argued in Ref. |^ that lattice calculations which 
linearly extrapolate the p mass as a function of mjr should 
receive negligible corrections from the two-pion interme- 
diate state. The calculations described here are in agree- 
ment with this conclusion, but as can be seen from the 
short-dashed curve in Fig. ^, the dependence of the real 
part of the vector self-energy in the region of the phys- 
ical pion mass is not smooth, showing the effects of the 
threshold crossed at 4m^. This rapid dependence on 
the mass of the intermediate states suggests that such 
extrapolations must be made with caution. 

The dashed curve in Fig. || shows the dependence of the 
self energy -m^p, m„, mi) when the mass of the in- 
termediate vector meson is my and the mass of the inter- 
mediate pseudoscalar meson mi is varied. Because of the 
different Lorentz-covariant structure of the vector- vector- 
pseudoscalar transition amplitude A^j^(pi,p2) given in 



TABLE II. Mass shifts to p and u) mesons arising from sev- 
eral two-meson intermediate states. Each channel provides 
less than a 10% correction to the total mass of the vector 
meson. The sum of these shifts, which is not observable, is 
given in the bottom row. The difference between these sums 
is m^j — rUp = 24.8 MeV, which can be compared to the ex- 
perimental value of 12 ± 1 MeV. 



channel 


Amp (MeV) 


Am^ (MeV) 


TTTT 


-69.8 




KK 


-14.8 


-14.8 


LUTV 


-22.5 




pvr 




-67.5 


ojr] 




-13.1 


PV 


-13.1 




K'K 


-11.9 


-11.9 


sum 


-132.1 


-107.3 



Eq. (2.8), and that the threshold for p — > lot: has not 
been reached even for mi = 0, a much weaker depen- 
dence on the mass mi is observed compared to that 
of n^^(— m^, mi, m2). Nonetheless, it may be suffi- 
ciently accurate to neglect pseudoscalar-vector interme- 
diate states with masses mi > vip when calculating the 
self energies, since these are smaller by at least a factor 
of two than the contribution from the lott intermediate 
state. 

It is interesting to note that the curves in Fig. || all ex- 
hibit an extremum for values of mi for which the sum of 
the intermediate-mesons masses (mi -I- m2) is within the 
range ^mp < (mi + m2) < mp. This is not coincidental, 
since nip is the value of the center-of-momentum energy 
at the decay threshold. For values of mi > nip — m2 the 
real part of the amplitude is a monotonic, decaying func- 
tion of mi. At threshold mi — nip — ni2, discontinuities 
in the second derivative of the real part of the self energy 
cause the rapid turnover observed in Fig. ||, so that the 
extremum of Ren(— mp, mi, m2) occurs for a value of mi 
slightly below the threshold value mi = nip — 7112- 

This suggests that in calculations of self energies for 
a particle of mass M , it may be sufficiently accurate to 
consider only those states for which the sum of the inter- 
mediate hadron masses mi -I- m2 lie on the range 
iM < m.i < M. This implies greatly simplified cal- 
culations of mixings between states where it is only neces- 
sary to consider the contributions from multiple-hadron 
states within this range of total masses. This is a general 
feature of such calculations and not a specific result of 
the particular channels considered herein, since it arises 
independently of the particular form and scales used to 
model the meson transition form factors f^^^ {pi,P2) 
and f^^P{pi,p2). 
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B. vector-meson mass shifts 

Results for the mass shifts due to the two-meson loops 
considered in this study are summarized in Table ||. It 
should be noted that only the difference between the total 
mass shifts of the p and u mesons is physically meaning- 
ful. The absolute size of each term does, however, give 
some indication of the importance of the dressing by two- 
meson loops. The largest mass shifts are that of the p 
mass due to the two-pion loop, and that of the lo due to 
the pTT loop, which are about —10% of the bare mass. 

Note that the p mass shift due to the two-pion loop 
is roughly —70 MeV, while the two-pion p mass shift in 
Ref. was between —10 and —20 MeV. The difference 
is due to an additional contribution to the vector self- 
energy which arises from a mrpp contact term. This term 
provides a contribution to the vector-meson self energy 
llfj,^{q) that is independent of q. In Ref. this term is 
added to provide current conservation for the vector p; 
i.e., that Il^^{q^ = 0) = 0. Since one can show that this 
term contributes equally to the self energies of the p and 
u mesons, it cannot contribute to their mass splitting 
and so has not been included here. 

To allow for the direct comparison with Ref. Q , the ef- 
fects of including this contact term into the present study 
can be reproduced by subtracting Il{q^ = 0) from the self 
energy n(g^). This procedure ensures that the condition 
obtained in Ref. ^ that Il^^{q'^ = 0) = is reproduced. 
The resulting m ass sh ift due to the two-pion loop ob- 
tained from Eq. (2.21) would then be approximately 



SuIt, 



Re 



{-^l ,m^,m^)-n^^(0,m^,TO^) 
2 mo 



« -9.2 MeV. (5.2) 
This mass shift for the p meson is similar to that of 

Ref. i. 

As pointed out in Ref. the strange- meson and 
vector-?7 intermediate states cannot contribute to the 
mass splitting; they are included here to illustrate that 
as the masses of the intermediate state mesons increase, 
the resulting shifts decrease rapidly. For example, the 
shift in the p mass due to the KK intermediate state is 
about one fifth of that due to two pions. Similarly, the 
shift in the uj mass due to the pr] intermediate state is 
roughly half the size of that due to tun. 

The result of adding the shifts due to the various inter- 
mediate states is that they largely cancel, as in Ref. 
to give a mass splitting of roughly 25 MeV, which is of the 
same order (on the scale of the bare mass) as the experi- 
mental value of 12 ± 1 MeV. Of course, a complete calcu- 
lation would have to consider all possible multiple-hadron 
intermediate states. It is encouraging that, within this 
model, a reasonable value of the mass splitting is achieved 
using only a handful of meson loops, and that the indi- 
vidual contributions due to additional loops are expected 
to be small compared to the dominant contributions eval- 
uated here. 




-0.6 -0.4 



FIG. 9. Ratios of the forward- forward (dashed curve), 
backward-backward (solid curve) and forward-backward or 
backward-forward (short-dashed long-dashed curve) time or- 
derings of the real part of the two-pion loop integral to 
the sum of all four, as a function of the vector-meson mass 
squared. 



C. time orderings 

Figure O shows the ratio of the forward- 
forward n[^+l(q^), backward-backward n[^^l(g^), and 
forward-backward II'^ '(9^) [equal to backward- forward 
nl ^1 (q-^)] time orderings of the real part of the two-pion 
loop integral to their sum n'^'^(g^), as a function of the 
vector-meson mass squared q^. As expected from the 
analysis in Sec. |l|, at = the forward-forward time 
ordering nt++l(g^ is equal to the backward-backward 
time ordering III '(9^)- It is also simple to show [by 
a change of variables from fc4 to —ki in Eq. |3.3| that 
n[+~l(g2) = n[-+I(g2) at aU values of g^. Although not 
shown here, an examination of the derivative of n[++l (g^) 
indicates a discontinuity at g^ = 4toJ indicating the on- 
set of a branch cut in 11'^'^ (g^) due to nt++l(g'^), which 
verifies the identification of this time ordering with the 
propagation of two forward-propagating pions. 

In Fig. H one observes that at g^ — — the forward- 
forward time ordering is dominant and the backward- 
backward contribution is only a few percent of the to- 
tal. This realizes the expected strong suppression of the 
backward-backward time ordering, due to the vertices 
and the large energy denominators for the ppmr interme- 
diate state. Note also that both the forward-forward and 
backward-backward contributions have the same sign, so 
that the two-pion loop contribution to the p self energy 
is not reduced by a cancellation between these two time 
orderings. This establishes that, in the present quantum 
field-theoretic framework, the dominant contribution to 
the p self energy is from two forward-propagating pions. 
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The form factor G^"" 
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FIG. 10. The p-meson electric-monopole form factor. The 
dashed curve is the resuh G|'((7^) from Ref. ^ that rises from 
quark-antiquark substructure of p meson. The solid curve is 
the real part of the form factor Ggf g ) that results when the 
effects of pion loops from Eq. ( 2.54 ) are added to the results 
from Rcf [i. 



and that the difference between the mass shifts obtained 
in this approach and those of the quantum-mechanical 
approach of Ref. Q is not due to the presence of addi- 
tional time orderings in a L orentz covariant framework. 
As discussed in Sec. 3.3, the forward-backward and 



backward-forward time orderings are not zero in the 
present framework; Fig. ^ shows that they also make neg- 
ligible contributions to 11'^'^ (q^) 



at = — TO^ 



D. vector- meson electromagnetic form factors 

The final application of vector-meson dressing by two- 
meson loops that is considered here is the correction to 
the EM f orm factors and charge radius of the p meson. In 
Sec. [IC, an expression for the pion contribution to the 
EM form factor G^'^(Q^) was obtained in terms of the 
same meson transition form factor f^^^{pi,P2), given 
in Eq. ( |4.9| ) with Eq. (4.11), that was used to obtain the 
TTTT and KK contributions to the p-meson self energy. 

In th e pre sent calculation of GJf (Q^), the final term 
in Eq. (2.54) (arising from the surface contribution of 
the integration by parts) is small and varies slowly with 
w 0, hence it is safe to neglect this term for the de- 
termination of the charge radius. The resulting G]^{Q'^) 
is well described for spacelike values of the photon mo- 
mentum squared < < 0.8 GeV^ by the form 



(5.3) 



with B « 1.39 and the value of A as determined from 



(Q^) is added to the quark- 
antiquark contribu tion Cj^jQ"^) obtained by Ref. ||] 
according to Eq. ( 2.55| ). The resulting form factor 
G^''+'^'^(Q2) is plotted as a solid curve in Fig. M It 
falls more rapidly with the photon momentum Q than 
does the form factor G^(Q^) (dashed curve) indicating 
that the inclusion of the tttt intermediate state into the 
p meson self energy has increased the charge radius of 
the p meson. The charge radius of the p meson due to 
its quark-antiquark substructure a nd its mixing with the 
two-pion state obtained using Eq. ( p^ ) IS 



0.67 fm. 



(5.4) 



This value may be compared to the value obtained in 
Ref. |, 



1 

V qql 



0.61 fm, 



(5.5) 



which results from only the quark-antiquark substructure 
of the p meson. One observes that including pion loops 
increases the charge radius of the p meson by approxi- 
mately 10%. This is consistent with the increase of less 
than 15% observed in a similar study of the pion-loop 
contribution to the tt meson charge radius 1^ . 

Although, the short lifetime of the p meson excludes 
the possibility of directly measuring the EM charge ra- 
dius of the p meson at present, one can define a "diffrac- 
tive radius" from the t dependence of diffractive p-meson 
electroproduction cross sections Q. As with the EM 
charge radius, the diffractive radius receives contribu- 
tions from its quark substructure as well as from mixings 
with multiple-hadron states. However, an important dis- 
tinction between the charge radius and the diffractive 
radius is that since diffraction arises from a strong inter- 
action, electromagnetically neutral particles, such as the 
Lo and p° mesons, will have diffractive radii even though 
they have no well-defined charge radius^. 

The results of this study suggest that the most signifi- 
cant contribution from two-meson loops to the diffractive 
radius of the p meson would be due to the two-pion loop. 
However, G-parity forbids the lo meson from mixing with 
the two pions, so that the diffractive radius of the uj me- 
son would be unaffected by the two-pion loop. Therefore, 
the diffractive radius of the w meson is expected to be 
smaller than that of the p meson. An observed differ- 
ence in the t-dependence of diffractive electroproduction 
of p and Lo mesons would be indicative of their differing 
diffractive radii, and such a measurement could provide 
a means to estimate the contributions to these radii from 
pion loops. 



Eq. (2.57) is A = SZ^^ = 0.134. 



* To be electroproduced diffractively, the vector meson must 
be able to mix with the photon, and so must be neutral (e.g., 
uj, p", 0, J/ip). Therefore, only for neutral vector mesons can 
one define a diffractive radius. 
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VI. CONCLUSIONS 

It is shown that, in a covariant model based on studies 
of the Schwinger-Dyson equation of QCD which assumes 
exact SU(3)-fiavor symmetry, the contributions to the 
self energies of the p and uj due to several pseudoscalar- 
pseudoscalar and pseudoscalar-vector meson loops are at 
most 10% of the bare mass. The result for the mass shift 
of the p meson due to the two-pion loop is in agreement 
with a previous calculation of this quantity using an effec- 
tive chiral Lagrangian approach. Such contributions are 
found to decrease rapidly as the mass of the intermediate 
mesons increases beyond mp/2. 

The mass shifts due to several two-meson intermediate 
states are compared with those from an extensive study 
within a nonrelativistic framework of the p-LO mass split- 
ting, and are found to be smaller, especially for interme- 
diate states involving higher-mass mesons. A mass split- 
ting of — nip Ri 25 MeV is found from the tttt, KK, 
ujir, p-K, LOT], prj and K* K channels. A complete calcula- 
tion which evaluates contributions from all possible mul- 
tiple hadron intermediate states and breaks SU(3)-flavor 
symmetry is beyond the scope of the present work. How- 
ever, these results suggest that such a calculation should 
exhibit rapid convergence as the number of two-meson 
intermediate states is increased by including states with 
higher masses. This implies that inclusion of two-meson 
loops into the vector-meson self-energy yields a small 
correction to the predominant valence quark-antiquark 
structure of the vector meson. 

The part of the vector-meson self energy which corre- 
sponds to two pions propagating forward in time has been 
shown to dominate a Lorentz-covariant calculation of the 
two-pion loop contribution. This implies that the addi- 
tional time-orderings of this loop, necessary to maintain 
Lorentz covariance, are not responsible for the reduced 
size of the mass shifts found in the present framework 
when compared to quantum-mechanical treatments. 

As a test of the self-consistency of this calculation, the 
contribution of the two-pion loop to the p meson EM 
form factor is evaluated, and is shown to provide a mod- 
est increase of about 10% to the charge radius of the p. 
Although this charge radius is difficult to observe, it may 
be possible to observe the effects of such an increase by 
examining the difference in the "diffractive radii" of the 
neutral p and u) found from the ^-dependence of diffrac- 
tive electroproduction cross sections. As the two-pion 
loop cannot contribute to the it may be possible to 
observe its effects on the p by measuring a difference in 
these diffractive radii. 
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